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ABSTRACT 

We present an analytical model which reproduces measured galaxy number counts from surveys in 
the wavelength range of 500 (im to 2 mm. The model involves a single high-rcdshift galaxy population 
with a Schcchtcr luminosity function which has been gravitationally lensed by galaxy clusters in the 
mass range 10 13 to 10 15 M©. This simple model reproduces both the low flux and the high flux end of 
the number counts reported by the BLAST, SCUBA, AzTEC and the SPT surveys. In particular, our 
model accounts for the most luminous galaxies detected by SPT as the result of high magnifications 
by galaxy clusters (magnification factors of 10-30). This interpretation implies that submillimcter and 
millimeter surveys of this population may prove to be a useful addition to ongoing cluster detection 
surveys. The model also implies that the bulk of submillimeter galaxies detected at wavelengths larger 
than 500 ism lie at redshifts greater than 2. 

Subject headings: galaxies: clusters — gravitational lensing — submillimcter: galaxies 



1. INTRODUCTION 

Over the last decade, submillimcter (submm) surveys 
have yielded significant advances in our understanding of 
the galaxy population responsible for the high-redshift 
component of the cosmic infrared background (CIB) 
dFixsen et al)[T996t|Smail et al|1997llHughes et ai1ft998f 
Barger et alJll998t iDwek et al l 11991 iFixsen et al.l | 1998t 
Greve et all 2004 iPope et al.l 120061: iCoppin et"aTll200a 
Devlin et all 120090 . With typical far-infrared (FIR) lu- 
minosities > 10 12 L©, submm galaxies are presumed to 
be the high-redshift counterparts to (ultra) luminous 
infrared galaxies (LIRGs, ULIRGs). The high lumi- 
nosity of these galaxies is the result of star formation 
rates of 100-1000 M Q yr _1 . Approximat ely half of these 
galax i es are located at 1.9 < z < 2.9 (jChapman et al.l 
120051 : lAretxaga et al.l 120071 ), dominating the total star 
formation rate at this ep och (jPerez-Gonzalez et al.ll20"o"5t 
iMicha lowski et al.l 



One way to express the results of submm surveys is 
through number counts of galaxies as a function of flux 
for each observed wavelength. The shape of these counts 
has been interpreted as arising from different popula- 
tions of g alaxies whose characteristic s evolve over cos- 
mic t ime (lLagache et ail 120031 l2004t iPearson fc Khanl 
120091: iLe Borgne et al.l 12009( 1 These empirical models 
have successfully reproduced the counts. However, they 
may be masking a simpler explanation for the depar- 
ture of the counts from a Schechter distribution at the 
high- flux end: magnific ation due to high redshift galaxy 
clusters and groups fe.g. lBi^ll99llPerrotta et al.1120021: 
INegrello et atH^CXTTh . 

Millimeter wavelength surveys have also aimed at de- 
tecting ga]axy_dujrters_vhi the Sunyaev Zel'dovich (SZ) 
effect (jHincks et al.l 120081: ICarlstrom et all I2009H : the 
first results, including CMB power s pectra and cluster 
catalogs, have been released recently dFowler et aLll2010l : 
iStaniszewski et all 120091 : iVanderlinde et al.l 12010( 1. The 
number of detected clusters remains relatively low, pri- 
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marily due to the low value of ag. However, other ef- 
fects could be reducing the sensitivity of the surveys (e.g. 
ILima et alll2009l l20loh . 

The South Pole Telescope (SPT) has measured num- 
ber counts of dusty galaxies at wavel engths A = 1.4 mm 
and 2.0 mm over an area of 87 deg 2 (jVieira et alll2009H . 
The observed numbers at the bright end are higher than 
expected: these galaxies are either at high redshifts and 
intrinsically exceptionally luminous, or have been magni- 
fied by gravitational lensing, or are simply at much lower 
redshifts than the bulk of the population of submm galax- 
ies. The latter possibility is disfavored by the lack of de- 
tected counterparts i n other surveys th at probed the low 
redshift population (jVieira et al.l l2009h . The possibility 
that these galaxies are at high redshifts and intrinsically 
bright would require them to be far more luminous than 
an underlying Schechter- like luminosity function would 
permit. Thus the favored explanation is that they have 
typical luminosities for high-z galaxies, but have been 
magnified by foreground galaxies or clusters. In fact, 
lensing of high-redshift background su bmm galaxies has 
been observed in a number of systems dSmail et allll997l 
20021: [Wilson et al.ll2008t iRex et al.ll2009HGonzalez et al l 
2009tlSwinbaiik et al.ll2010f ). 



In this Letter, we explore the possibility that the exist- 
ing observed galaxy number counts over a wide range of 
wavelengths can be reproduced by a single population of 
galaxies at high-redshift. Foreground galaxy groups and 
clusters gr avitationally lense s the background submm 
population (jLima et alll2009l ) and leads to significant en- 
hancements of the high-flux end to the galaxy counts. 
In § [2] we describe the lensing magnification formalism, 
which we then apply to a high-z galaxy population and 
present results in §[3] We discuss implications for high-z 
galaxies and the cluster searches in § |4j 

Throughout, we use a fiducial cosmology for a flat uni- 
verse with parameter values based on the results of the 
Wilkinson Microw ave Anisotropy Probe fifth year data 
release (WMAP5. Ikomatsu et all 120091 ) . The cosmo- 
logical parameters (and their values) are the normaliza- 
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Figure 1. Intrinsic and lensed number counts dn/dS for a 
Schechter function describing galaxies at different redshifts. Also 
shown are the observed counts for SPT dusty submm galaxies at 
A = 2.0 mm, after removal of low-redshift galaxies with IRAS coun- 
terparts. 

tion of the initial curvature spectrum 8q(= 2.41 x 10~ 4 ) 
at k = 0.02 Mpc -1 (corresponding to as = 0.8), its 
tilt n(= 0.96), the baryon density relative to critical 
Q b h 2 {= 0.023), the matter density n m h 2 (= 0.13), and 
two dark energy parameters: its density f2rjE(= 0.74) 
and equation of state w(= —1), which we assume to be 
constant. Since lensing effects depend on the abundance 
of dark matter halos, which is exponentially sensitive to 
linear matter perturbations, we also consider changes in 
as consistent with the WMAP5 errors of Acts w 0.03. 
Our central value and unc ertainty for cr 8 is con sistent 
with the WMAP7 analysis (jKomatsu et al.ll2010H . 

2. NUMBER COUNTS WITH LENSING MAGNIFICATION 

In recent papers (jLima et al.ll2009t Uain fc Limall2010D 

we have presented a halo model for calculating the effects 
of lensing magnification by galaxy groups and clusters. 
Here we specialize to the case of steep galaxy counts at 
high redshifts, where lensing eff ects are quite dr amatic. 
We assume a Schechter function (|Schechterlll976|) for the 
intrinsic number density distribution of a population of 
galaxies 



!tt = !?L(JL) e s/s- 



dS S* \S* 



(1) 



where n* , S* and a are free parameters. Due to lensing 
magnification by intervening halos, the intrinsic dn/dS 
is changed to its observed counterpart as 



dn bs(S bs) 

dS b: 



= / dfj, 



P(n) dn ( Sobs 



/i dS \ /i 



(2) 



where fi is the lensing magnification and P(/i) is its 
probability for a given galaxy population at rcdshift z s . 
Conditional probabilities quantify the effects of different 



magnification ranges on the observed flux density Sobs- 
The integrand of Eq. [2]dcfhics the probability P(^|S Q bs) 



( rfn bs(<S , obs)\ 1 P(aO dn ( S, 



\ dS bs 



// dS \ n 



Jobs 



,(3) 



which can be interpreted as the relative co ntribution of 
a given \i to the total G?n Q bs / dS b s at Sobs (jPaciga et al.l 
USSH). Similarly P(/x min |S bs) = /^. n dfi P(/i|S bs) mea- 
sures the integrated contribution from all fi > /i m m- The 
mean magnification at a given Sobs is defined as 



(ft) (Sobs) 



d(i [i P(/i|S bs) 



(4) 



The distribution P(fJ.) can be esti mated either by ra 
tracing on N-body simulations (e.g. IHilbert et aLll200 
or by semi-analy tical methods (e.g. iPerrotta et al.ll2002l 
iLima et al.ll2009h . integrating halo contributions on the 
line of sight up to the source redshift 



P(> M ) 



dz, 



D\(zi) 
H(zi) 



A/th 



d]nM dn(z h M) 
d\nM 

(5) 

where Da is the angular diameter distance, H is the 
Hubble parameter, dn/dlaM is the halo mass- function 
and Ail^ = A^l fl (z s , zi, M) is the cross-section for mag- 
nifications larger than \i produced by halos of mass M 
at redshift z\ on sources at redshift z s . The integrand 



d?P(> fi) _ D\(zi) dn(z h M) 
d\nM dzi ~ H(zi) d\nM 



An, 



(6) 



gives the range of halo masses and redshifts contributing 
most to the probability of a specified minimum magnifi- 
cation fj.. 

In summary, Eqs. [5] and [S] give the total effect on the 
counts, Eq. [3] indicates which magnifications contribute 
most to the given S bs and Eq. [6] tells us which halo 
masses and redshifts contribute to a given magnification. 

We use this halo-model P(/x) and c orrect it for anum- 
ber of effects. First, as described in ILima et al.l ()2009f ) . 
we match our P(/x) at large magnific ations to that of ray - 
tracing in dark matter simulations (|Hilbert et al.|[2007|) 
by tuning the ellipticity of our halos. Next, we account 
for the effect of luminous matter, which can lead to 
higher densities via gas coo ling, using the simulation re- 
sults of IHilbert et al.1 (j2008l ) . We also correct for the com- 
bination of finite source size and multiple image effects. 
Magnification effects, especially for large magnifications, 
are sensitive to the value of as since it affects the abun- 
dance of cluster halos. In the next section we discuss how 
we account for the uncertainties in our model by giving 
a range for our predictions. 

Our analytical calculation of P(n) has some advan- 
tages over the approach of numerical simulations (we can 
easily study changes in source redshift, as and the con- 
tribution from different halo masses and redshift), but it 
also has some limitations. We only use the one-halo term, 
which is accurate at the high magnifications relevant for 
the effects considered here but overestimates the lensing 
contribution at fi ~ 1. We do not include a distribu- 
tion of ellipticities or halo substructure, which can also 
increase magnification cross-sections. We have instead 
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Figure 2. (Left) : Intrinsic and lensed dn/dS for a Schcchtcr function describing galaxies at z s = 3.0 and different wavelengths. Dashed 
lines indicate the intrinsic Schechter functions (with different S* values) and the dark shaded regions display the range of lensing predictions, 
as described in the text. Also shown are observed counts for BLAST at A = 500 /^m, SCUBA at A = 850 /^m, AzTEC at A = 1.1 mm 
and SPT at both A = 1.4 mm and 2.0 mm. The SPT counts are for dusty galaxies, after removal of galaxies with synchrotron emission 
and galaxies with IRAS counterparts. No similar removal has been applied to the BLAST data, which includes both high and low redshift 
galaxies. Notice that we do not display the lensing predictions for SCUBA and AzTEC, since these data do not require lensing in the 
measured fluxes. A prediction for MUSTANG at A = 3.3 mm is also shown in the light shaded region. (Right) : Unified scaled curves 
show'mgdh/ dS and the various data points. The only parameter used in the scaling is S 1 *; its values at the different wavelengths are shown 
in Fig. [3] 

tuned the average halo ellipticity to match the P(n) mea - 
sure d in dark m a tter s imulations of lHilbert et al.l (|2007f ) 
(see iLima et al.1 (|2009fl for a detailed discussion). And 
whereas we account for t he effects of ba r yons observed in 
numerical simulations bv lHilbert et al.1 < 2008T ) . these au- 
thors note that their simulation still underestimate bary- 
onic effects for halos of smaller masses as they do not 
predict sufficient numbers of multiply imaged quasars. 
Finally, the effect of finite source size is very uncertain 
given the lack of our knowledge about submm galaxies 
and its sens itivity to the pr ecise caustic structure of the 
lenses (e.g. iLi et al.l (|2005T i ). A proper inclusion of all 
these missing effects would likely increase the magnifica- 
tion probabilities compared to our current model. 
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3. RESULTS 

In Fig. [1] we illustrate the lensing effect on an intrinsic 
Schcchtcr function distribution given by Eq.[T]for sources 
at different redshifts. In all results presented here, we 
fix a — —1.0 and n* = 5 x 10 3 dcg~ 2 . Changing to 
a = —1.5 does not have a significant effect - while it 
m atches the faint end b ehavior of dn/dS for the model 
of lLagache et al.1 (|2004f ) . we preferred to use a = —1.0 
since this fits better the number counts at shorter wave- 
lengths. With these pa rameter values , our c ounts are 
lower than the model of lLagache et al.l (|2004l ) at all S, 
which also ensures that the total flux does not exceed the 
CIB (jDwek et al.lll998L iFixsen et allfloM ). 

Fig. [5] illustrates our main results: we show pre- 
dicted number counts that include lensing (gray bands), 
assuming galaxies at z s = 3.0, along with measured 
number counts from submm surveys at different wave- 
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Figure 3. The flux S* in the Schechter function is shown for dif- 
ferent surveys at various wavelengths (symbols). The three curves 
show the expected scaling for submm galaxy SEDs at z s = 3.0 for 
different values of spectral index /3 and temperature T. The curves 
are normalized at the value of S* for A = 850/xm. 

lengths. As indicated in the panel s , thes e are: BLAST 
at A = 500 urn (iDevlin et al l [2P09), SCUBA at 
A = 850 fim (jCoppin et al.l I2006D . AzTEC at 1.1 mm 
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Figure 4. Probability P(/J-mln\Soba) of the minimum magnification fi m i n , given an observed flux density 5 b s (left panel), and the 
corresponding average magnification (fi) (right panel). We assume a Schechter function describing galaxies at z s = 3.0 which, after lensing, 
predicts counts consistent with those of SPT dusty submm galaxies at A = 1.4 mm (see Fig. [2j. 



(jAustermann et al~ll2010( ) and SPT dusty submm galax- 
ies at A = 1.4 mm and 2.0 mm dVieira et al.|[2009h . In all 
ou r results, the SPT number counts correspond to those 
of iVieira et all (|2009l ). after removing both synchrotron 
emission galaxies as well as low-redshift galaxies that 
have matches with gal axies in the Infrared Astronomy 
Satel lite survey (IRAS IMoshir et ail 119921: iFisher et all 
Ill^l01iverltal1l9m Predictions are shown by bands 
rather than curves to reflect the uncertainties in the 
model as discussed below. The intrinsic luminosity func- 
tion is the Schechter function described above. All the 
data sets can be fit by changing the single parameter S* 
once lensing magnification is included. This remarkable 
result implies that, within the measurement and theoret- 
ical uncertainties, a single high-z population of galaxies 
is sufficient to describe all the observations. The high 
flux measurements of BLAST and SPT are fit by highly 
magnified galaxies - if these counts were dominated by 
a population of a different galaxy type, it would be a 
coincidence that their relative counts fit the same scal- 
ing with wavelength as the fainter ("normal") popula- 
tion. Finally, we also show the prediction for the Mul- 
tiplexed Sojjid_TT^A I ray at 90 GHz (3.3mm) (MUS- 
TANG, iMason et all 120061 ) which has begun operating 
on the Green Bank Telescope. 

The lower bound for the predictions in Fig. [2] uses 
o~ g = 0.77 while the upper bound uses cr 8 = 0.83 - 
these reflect the uncertainties in the WMAP results as 
discussed above. The enhancement due to baryons is 
factor ed into these predictions following iHilbert et al.1 
(|2008| ) , though it is likely to be an underestimate of bary- 
onic effects as discussed above in Section 2. The effect of 
source sizes and multiple imaging is uncertain; we sim- 
ply assume that due to the finite source size the effective 
magnification is reduced by 50% to 25% (for the lower 
and upper bounds respectively) 



The right panel of Fig. [2] shows all points rescaled by 
plotting 



where n = n/n* and S = S/S*. In Fig. [3] we show 
values of S* used for each wavelength. We compared the 
frequency scaling of S* with that of a typical Spectral 
Energy Distribution (SED) of submm galaxies SED(A) cx 
e(X)B(T, A), redshifted to z s — 3 with emissivity e(A) = 
1 — exp[(— Ao/A)' 3 ] and blackbody spectrum B(T,X) at 
temperature T. The values of S* arc consistent with 
,8 = 1 — 2 and T = 30 — 40 K, as expected for high 
rcdshift submm galaxies. 

The fit to the BLAST counts at A = 500, 350 and 
250 yum falls further below the high flux measurement at 
shorter wavelengths, sugg esting the need for a lower red- 
shift population. Indeed, lEales et al.l (|2009f ) have iden- 
tified the radio and 24 (im counterparts of the bright 
BLAST sources. At 250 (im almost all of the bright 
sources are identified as being at redshifts lower than 1. 
One third to one half of the sources in the highest flux 
bin at 500 /zm come from sources with z < 1. Remov- 
ing these would lower the corresponding point in Fig. [3J 
in agreement with the model. The remaining sources at 
500 fj,m (and less than a tenth of the sources at 250 /im) 
are likely the result of lensing. Similar results should 
be expected with the upcoming release of the large-area 
Herschel surveys. 

Since the lensed distributions at z s = 3.0 arc con- 
sistent with SPT data points at both wavelengths, we 
study the range of magnifications and halo masses that 
would be contributing most in this case. As we consider 
larger values of Sobs, in particular S oos /S* ^ 10, the ob- 
served sources come from intrinsically low flux sources 
which have been magnified significantly. Fig. [4] shows 
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Figure 5. The integrand d 2 P(> fi) /d In Mdzi as a function of 
halo mass M for different values of fi and Z[. We have arbitrarily 
normalized these curves for better visualization. Note that the 
peak lensing contribution for source galaxies at z s = 3 comes from 
lens halos at Z\ ~ 0.5. This does not include baryonic effects, which 
boost the contribution from M ^ 10 13 h~ 1 Mq halos. 

-P^minlSobs) and (fx) as a function of Sobs for SPT (A = 
1.4 mm) and indicates magnifications that contribute 
most at each S bs- For instance, for S bs = 20 — 40 mJy, 
the right panel shows that (/i) ~ 20 — 30. Another way 
to see this is via -P(/i m i n |S bs) (1°^ panel), which inte- 
grates out the effects above a certain fi m i n , and shows 
that for S'obs = 20 — 40 mJy, we have P(Li m in\S bs) > 0-5 
for fi m i n ~ 10 — 20. 

These results imply that magnifications of 10 — 30 
are necessary to explain the boost in dn/dS at S'obs ~ 
10 — 40 mJy, if it is due to lensing of an intrinsic 
Schechter distribution. Note that due to the finite size of 
submm galaxies, their magnifications must have a cut-off, 
which has been estim ated to be in the range ll ~ 10 — 40 
(jPerrotta et al.l[2002l ) for galaxy lenses, and is probably a 
factor of two or so larger for more massive lenses. Indeed 
galaxies have been mea sured with estimated magnifica- 
tions of at least ~ 45 dPaciga et al.l 120081 iKneib et al.l 
12001 . 

In Fig. [5] we show d 2 P/d\nMdzi as a function of halo 
mass for different values of /i m ; n and z\. This indicates 
that halo masses above 1O 13 /i _1 M0 contribute signifi- 
cantly to magnifications of 10 — 30, with most of the 
contribution coming from ~ 10 14 /i _1 M Q . Note that this 
does not include baryonic effects, which boost the con- 
tribution from lower mass ( iS 10 13 Ii~ 1 Mq) halos as dis- 
cussed in Section 2. 

We have assumed that all source galaxies are at a fixed 
redshift z s , whereas in reality they have a redshift dis- 
tribution which needs to be incorporated in the compu- 
tation. The curves at different source rcdshifts shown in 
Fig. [T] provide approximate limits for what we can ex- 
pect from a redshift distribution of source galaxies. Our 
results imply that the bulk of the population of submil- 



limeter galaxies are at rcdshifts z £ 2. The alternative 
explanation for the high flux measurements is to have a 
significant fraction of galaxies at very low z, which would 
easily have been observed in surveys such as IRAS (in- 
deed the SPT data shown remove the small fraction of 
such sources) . And as shown in Fig. [TJ galaxies at higher 
redshift but still at z ;$ 1 would not match the measure- 
ments as the magnification boost is insufficient (again 
within the context of an underlying Schechter distribu- 
tion). 

4. DISCUSSION 

We have considered the possibility that bright millime- 
ter and submillimeter galaxy counts arise largely from a 
galaxy population at high rcdshifts that is lensed by in- 
tervening galaxy groups and clusters. Our model predic- 
tions match the counts from 500 /im to 2 mm from the 
BLAST, SCUBA, AzTEC and SPT surveys (sec Fig. [2). 
We find that the high flux SPT number counts can be 
explained by highly magnified galaxies from this high- 
z population (once the known, low-z counterparts de- 
tected in IRAS are removed). This high-z galaxy popula- 
tion is described by a Schechter luminosity function with 
L* ~ 2.5 x 1O 12 L and a source redshift z s = 3.0. Our 
model predictions fit the data for 500 /im < A < 2 mm by 
varying S* with wavelength within the range of typical 
submillimeter galaxy SEDs (see Fig. [3]). 

Our model has some simplifying assumptions, such as 
the fixed source redshift z s = 3, and there are significant 
measurement and theoretical uncertainties. A complete 
analysis would require a more detailed treatment of the 
lensing and inclusion of the measured galaxy clustering. 
In addition, it has been established that at the shorter 
wavelengths probed by BLAST, an increasing fraction of 
sources lie at low-z - hence we can expect some smooth 
variation in the fraction of low-z galaxies with observed 
wavelength. 

Nevertheless our results in Fig. [2] imply that current 
number counts do not require an additional galaxy pop- 
ulation to explain the high flux measurements. Such a 
population has been invoked in theoretical models (e.g. 
lLagache et al.1 l20"03h and suggested as a possible expla- 
nation (as well as lensing of th e high-z populatio n) for 
the recent SPT measurements (|Vieira et al.l |2009() . In- 
deed, adding a significant fraction of a second popula- 
tion could cause the predictions to exceed the measured 
counts once magnification effects for the high-z popula- 
tion are included. It would also be difficult to explain 
how the scaling with wavelength of the high flux num- 
ber counts is the same as the lower flux counts if they 
came from different galaxy populations. Since lensing 
does not depend on frequency, our model naturally fol- 
lows this common scaling. 

Current SZ surveys with sensitivities of 3 — 7 x 
1O 14 /i _1 M cannot detect most halos that produce this 
lensing contribution. Conversely, looking for extremely 
bright objects in millimeter and submillimeter wave- 
lengths provides a way to find high-z lensing halos as- 
sociated with galaxy groups and clusters. Note that the 
number of halos that can be found in this way is only a 
small fraction of all halos within a given mass range. As 
discussed above, our model most likely underestimates 
the contribution from halos with M ;$ 10 13 M Q . It is 
therefore of great interest to investigate the lenses corre- 
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sponding to the bright sources in current data. 

Follow up observations with optical telescopes should 
be able to identify lensing group/cluster candidates up 
to z ~ 1. These clusters and groups host Brightest Clus- 
ter Galaxies with lumin osities L ~ 10 11 — 10 L w , based 
on the low-z results of Uohnston et all ()2007l ) from the 
Sloan Digital Sky Survey. Multi-band optical imaging 
with limiting magnitude of 23 — 24 (for the r band) would 
enable identification of these groups and clusters. Con- 
versely, targeted observations of the high magnification 
regions of known strong lensing clusters could provide 
detections of faint submm galaxies (which would lie be- 
low the detection threshold without the magnification 
boost). Similar to the use of clusters as gravitational 
telescopes in optical imaging, this may also help resolve 
submm galaxies. Planned observations with AzTEC and 
the Large Millimeter Telescope have considered such an 
approach (David Hughes, private communication). 

We thank David Hughes, Roxana Lupu, Joaquin 
Vieira, Kim Scott, Ian Smail, Eric Switzer and various 
members of the ACT collaboration for useful discussions. 
We benefited from discussions on strong lensing effects 
with Matthias Bartelmann, Gary Bernstein, Neal Dalai 
and Ravi Shcth. We are very grateful to Stefan Hilbert 
for sharing his simulation results and Eric Switzer for 
providing the SPT data. This work was supported in 
part by an NSF-PIRE grant and AST-0607667. 
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